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Abstract. Geometric phases of scattering states in a ring geometry are studied based 
on a variant of the adiabatic theorem. Three time scales, i.e., the adiabatic period, 
the system time and the dwell time, associated with adiabatic scattering in a ring 
geometry plays a crucial role in determining geometric phases, in contrast to only two 
time scales, i.e., the adiabatic period and the dwell time, in an open system. We derive 
a formula connecting the gauge invariant geometric phases acquired by time-reversed 
scattering states and the circulating (pumping) current. A numerical calculation shows 
that the effect of the geometric phases is observable in a nanoscale electronic device. 
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1. Introduction 

The study of geometric phases continues to be an intriguing subject. Since Berry's 
original work [1], various generalizations have been proposed [2, 3] to adapt to different 
physical applications in diverse fields. In Ref. [4] the concept of geometric phases for 
scattering states was introduced which allows to give a geometric interpretation of the 
formalism developed by Brouwer [5] for quantum adiabatic pumping in open systems. 
Specifically, Brouwer presented a compact formula for the pumped charge (current) in 
terms of the parametric derivatives of the time-dependent scattering matrix subjected to 
the oscillating potential, which has been identified as the geometric phase accompanying 
the scattering state associated with the time-reversed Hamiltonian. Quantum adiabatic 
pumping in open systems is subject to intense study [5, 6, 4, 7] due to the observation 
of adiabatic electron transport through a quantum dot subject to slow cyclic variation 
of gate voltages [8] . 

An intriguing question has been raised recently about the possibility to generate 
a circulating (pumping) current in a ring geometry by adiabatically varying external 
parameters [9, 10]. In Ref. [9], Cohen considered open systems as some subtle limit 
of closed systems and demonstrated that it is possible to reproduce the formulas by 
Landauer [11] and Buttiker, Pretre and Thomas [12] from the Kubo formula whereas 
Moskalets and Buttiker [10] examined this problem in the slow frequency limit of the 
Floquet scattering theory. These differences in approaches mean that the results are 
different, although they aimed at describing the same system. Such an unsatisfying 
situation calls for a thorough analysis of the problem. In this Letter, we address 
this issue from gauge field and geometric phase perspectives. We derive a compact 
formula connecting the circulating (pumping) current [13] and the gauge invariant 
geometric phases acquired by time-reversed scattering states. This is achieved by a 
purely geometric argument, based on the fact that pumped charge is additive for two 
consecutive pumping cycles. A generic feature of pumping in a ring geometry is that the 
momentum is time-dependent and so are the discrete energy levels. The closed system is 
different from the open in another aspect. In an open system there are only two different 
time scales, whereas in the closed ring geometry three different time scales appear, the 
adiabatic period, the system time and the dwell time. The physics behind this is the 
following. In the adiabatic limit, electrons are in an instantaneous scattering state. A 
small deformation in the scattering potential results in a small increment of geometric 
phases, which electrons can pick up at the dwell time scale. The dwell time being the 
time it takes for an electron to complete one scattering event. Since the electrons travel 
in a ring geometry, they have to obey the periodic boundary condition to which they 
adjust at the system time scale, t s . All electronic wave functions are extended over the 
entire ring thus inducing a circulating current. The latter is experimentally measurable 
in terms of a SQUID that measures the magnetic field produced by the circulating 
current. 
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2. The system 

Consider a quantum dot with two leads attached such that the two single-channel 
leads are bent back to form a ring [see Fig. 1(a)]. Suppose the dot is characterized 




(a) (b) 



Figure 1. (a) shows the setup with a time dependent scattercr, defined by its 
scattering matrix, S(t). Incident, A, and reflected wave, B. (b) shows the time 
scales relevant to the problem. Scattering takes place on r s , the dwell time (the time 
it takes for the scattering event) is t^. t s designates the time passing for the scattering 
event, it is only valid on the r s -scale. The time it takes for the wave function to adjust 
the boundary conditions to the scatterer, of the order of r s , must be shorter than the 
adiabatic changes in the scatterer, happening on the scale r. 

by a scattering matrix S(t) which depends on time t via a set of independent external 
parameters X = (X 1 , ■ ■ ■ , X u , ■ ■ ■ , X p ) oscillating slowly with frequency u. For a spin- 
independent scatterer, S(t) is a 2 x 2 matrix. Define the row vectors n L = (r,t), and 
n R = (t',r'), the unitarity of the scattering matrix implies that they are orthonormal 
n « . n /3 _ where a,f3 — L,R. From the boundary condition of the matching of 
the wavefunction for a ring geometry, the dispersion equation is written in terms of 
scattering matrices, det(SS w — 1) = with S w = exp(ikL)a x , where a x is the Pauli 
matrix. S w is the scattering matrix describing the remaining part of the ring except 
the dot. The energy spectrum is then discrete in the ring geometry. The adiabaticity 
requires that the frequency u> is small enough compared with the (averaged) level spacing 
A. On the other hand, for the instantaneous scattering matrix to make sense, it is 
necessary that the system time r s = L/v (v- electron velocity) is much greater than the 
dwell time Td[14] during which electrons scatter off the scatterer. That is, in the ring 
geometry, there are three different time scales: the period r = 2tt/uj, the system time 
r s = hA^ 1 , and the dwell time, r^. In the scattering approach to adiabatic quantum 
pumping in the ring geometry, we assume r >> r s » This condition ensures that 
it is legitimate to speak of instantaneous scattering states, consistent with Heisenberg 
uncertainty principle [15]. We illustrate this in Fig. 1(b) and notice that it is different 
from the adiabaticity condition discussed in Ref. [10]. Generically, the allowed (discrete) 
fc's, which in turn determine the discrete energy levels, are time-dependent due to the 
fact that the external parameters defining the scatterer X = (X 1 , • • • , X v , • • • , X p ) are 
changing with time and the scattering matrix are decisive for the dispersion, a feature 
different from scattering in open systems. 
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Discretized approach and the continuous limit. Suppose the system undergoes an 
adiabatic cycle with the period r, which is characterized by a closed loop in parameter 
space. Suppose that initially an electron is in a given discrete energy level, the 
adiabaticity ensures that it remains in the same level during the entire period. To 
take into account the time scales properly, we adopt a discretized approach, i.e., the 
entire interval [0, r] is divided into N pieces. Correspondingly, the closed loop may be 
regarded approximately as a polygon [see Fig. 2(a)]. Here we stress that each subinterval 




Figure 2. (a) displays the discretization of the adiabatic loop. Doing this, we can 
ensure that the scattering matrix formalism gives a correct description and distinguish 
the time scales associated with the theory. Here, the entire interval [0, r] is divided 
into N pieces. Each subinterval is of the same order of magnitude as the 

system time t s . At each instant, ij, two new vectors nf and nf are defined in terms 
of the instantaneous 5-matrix. (1)) displays the property that the geometric; phase 
is additive, and the pumped charge being a function of the geometric phase has to 
be additive. The lower contour displays a special case in which two cycles share the 
same initial configuration, represented by a dot, whereas the upper contour shows two 
separate cycles joined by two adjacent but opposite lines. 

[tj,tj + i] is at the same order of magnitude as the system time r s . At each instant U, we 
have two row vectors nf and nf . Since the argument is parallel, we drop the superscripts 
L and R for brevity. For two consecutive instants tj and , we may define the relative 
phase between the two vectors n« and iij+i by exp(i0j) = n^* • nj + i/|nj • n* +1 |. That is, 
4>i = Imlnnj* ■ n i+1 . Therefore we have the total phase for the closed loop (polygon) 

N N 

(f> = ^<Pi = ImlnJJni* ■ n m . (1) 

i=l i=l 
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That is, both n L and n R acquire a geometric phase <p L and <p R , respective, if the system 
returns to the initial configuration after one adiabatic cycle. Indeed, Eq. (1) is invariant 
under the (local) gauge transformation = e* ai rij, which arises from the fact that the 
absolute phase is not observable in quantum mechanics. In the continuum limit, we 
have 



This amounts to the statement that the gauge potential A is Im(n* • dn) up to an 
(integrable) immaterial term. Indeed, in the adiabatic limit, it is reasonable to assume 
that the gauge potential A only depends on n, dn, k, dk and L. From this we see that the 
most general gauge potential we can construct takes the form A = n* • dn + g(kL)dkL 
with g an arbitrary function, if we require that A transforms as 



under the gauge transformation. Here the second term proportional to dkL comes from 
the dependence of discrete fc's on time. Below we will see that g = i/4. Although the 
gauge potential A is formally identical to its counterpart in the case of open systems 
when expressed in terms of the row vector n, the explicit dependence of n on k which 
varies with time characterizes the difference. 

Similarly, we may define two column vectors h L = (r, t'), and h R = (t,r'), which 
are orthonormal due to the unitarity of the S matrix. Repeating the above argument, 
we conclude that n L and acquire, respectively, geometric phases <p L and (f> R , which 
take the form, = Im § n* • dn. Here the hat denotes the time-reversal operation since 
the row vectors and column vectors are connected via time reversal operation. 

Consistency with the periodic boundary condition. The above discussion indicates 
that there is a gauge group Ul{1) x Ur(1), which is necessary to accommodate the 
geometric phases 4>l and 4>r acquired by the two row vectors n L and iir, respectively. 
As we see, this results from the gauge freedom associated with the incident waves. 
In contrast to adiabatic pumping in open systems, one has to keep in mind that 
the scattering wave functions should satisfy the periodic boundary condition in the 
ring geometry. If we set ^ l{ x ) — Aexp(ikx) + B exp(-ikx) for —L/2 < x < 
and ^r(x) = C(ikx) + Dexp(-ikx) for < x < L/2, then the incident waves 
and the scattered waves are connected via B = Ar + Dt' and C = At + Dr' . On 
the other hand, the period boundary condition requires that B = D exp(-ikL), and 
C = Aexp(-ikL). Besides the dispersion equation which results in discrete energy 
spectrum, we have D = Ar / (exp(-ikL) — t'). Formally, we write the scattering state 
■0 as ij} — Aipi + DipR with ipi = exp(ikx) + r exp(-ikx) + texp(ikx) and ipR = 
exp(-ikx) + r' exp(ikx) + t'exp(— ikx). Then the Ul(\) x Ur{1) gauge transformation 
takes the form i/)l(r) ~^ ex P(^ Q; L(i : ?))V ; L(-R) ; which induces transformations on A and 
D: A — > exp(— %oli)A and D — > exp(— ian)D. Then we see that under such a gauge 
transformation the dispersion equation remains the same, but an extra phase factor 
exp(— i(otL — or)) appears in the relation between A and D. With tpi and ipR as the 




(2) 



A' = 



ida + A, 
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local basis, the time evolution governed by the Schrodinger equation induces parallel 
transport [16]. Taking into account the causality condition which states that scattered 
waves appear only after incident waves hit the scatterer, we may derive the expression 
of the gauge potential A (where g is determined to be i/4), resulting in the geometric 
phases 4>l and 4>r. When the system is brought back to the initial configuration, the 
scattering states ipL and i/jr acquire, respectively, geometric phases 4>l and 4>r- Now we 
may picturize the observable effect of geometric phases as follows. Suppose the system 
undergoes an adiabatic cycle. During this process, at each instant, electrons are in an 
instantaneous scattering state. When electrons feel a slow deformation in the scattering 
potential, they pick up a small increment of geometric phases at the dwell time scale, 
and then take time at the system time r s scale to redistribute themselves to maintain 
the periodic boundary conditions, thus inducing the circulating current in the pumping 
experiment setup. 

In the above we focused on scattering geometric phases acquired by scattering 
states. However, in the pumping setup, physical observables are connected with pumping 
geometric phases associated with time-reversed scattering states, a fact already known 
for pumping in open systems. In fact, there is another gauge group Z7x,(l) x Ur(1) 
associated with the gauge freedom of the scattered waves, which accommodates the 
geometric phases 4>l and (f)R. 

Connection between physical observables and geometric phases. Now let us establish 
the connection between the geometric phases 4>l and 4>r and the physical observables 
which are the charge Q a accumulated inside the dot and the charge Q c circulating in 
the ring during the entire pumping period r for a given occupied level. Suppose Q L 
and Q R are, respective, the charges pumped from the dot to the left and right sides. 
Obviously, we have Q a = — (Q L + Q R ). As for Q c , one may expect that it is proportional 
to (Q L — Q R ), i.e., . Q c = i{Q L — Q R ) with £ some undetermined constant. To determine 
£, we consider the special case in which no charge is accumulated inside the dot. In such 
a case, we have Q L = —Q R = Q. That is, the amount of charge Q is pumped from the 
left side to the right side of the dot. Then we have £ = 1/2. 

As observables, Q L and Q R must be gauge invariant, so they must be some functions 
of the geometric phases <j) L and (j) R . Since the left and right sides are independent, we 
have Q = f{4>) for a given closed loop (here we drop the superscripts for brevity). 
To determine the function /, we notice that the charge is additive. That is, if we 
consider a closed loop consisting of two consecutive closed loops C\ and C 2 [see the 
lower contour in Fig. 2(b)], then we have Q c = Qci + Qc 2 - On the other hand, the 
geometric phase is Abelian, implying that <pc = <f>Ci + 4 > c 2 - That is, /(0Ci + 0c 2 ) — 
fi&Ci) + /(0c 2 )- /.From this we see that the function must be linear, i.e., Q = c<j) 
with c an undetermined constant. Combining this with the Friedel sum rule which 
states that charge accumulated inside the scatterer follows SQ a = e/(27ri)<51n(detS'), we 
have c = — l/(27r). Actually, the argument is applicable to any two arbitrary loops as 
illustrated in the upper contour in Fig. 2(b). 
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In the context of adiabatic pumping, we are interested in pumping cycles without 
any charge accumulated inside the dot. In this case, we have 4> L = —<j) R . So we may 
define the circulating current I as I = Q c /t, which takes the form 

(3) 

47TT 

/,From this we immediately conclude that for an embedded quantum dot with the mirror 
symmetry, i.e., the left-right symmetry, there is no pumping circulating current. This 
is consistent with Moskalets and Biittiker [10] who stated that the spatial asymmetry 
of the scatterer is a necessary condition for the existence of an adiabatic pump effect. 

We emphasize that Eq. (3) describes the contribution of a certain discrete energy 
level to the circulating current. The total current results from summing up all occupied 
energy levels, which depends on the number of electrons in the ring. 

A numerical example. Consider a quantum dot modeled by a potential V(x), which 
is defined as for \x\ > a, V\ for —a < x < —b, V 2 for \ x \ < b, and V3 for b < x < a. Here 
x denotes the coordinate along the ring. The same potential has been used to model the 
quantum dot embedded in a double path interferometer proposed to directly observe 
scattering geometric phase [7]. Then the instantaneous spin-independent 2x2 scattering 
matrix S(t) for the dot may be determined from the solution of the Schrodinger equation 
(—(h 2 /2m)d 2 /dx 2 + V(x) — E)ip = 0. For a dot of size 800 nm the energy level spacing 
is of the order of 4.5 meV. The Coulomb energy, assuming a dielectric constant of 10, 
is of the order of 0.08 meV. Thus, the dimension of the dot is such that the Coulomb 
energy is much less than the separation between the resonances and can be ignored. 
Also the spin-dependent scattering inside the dot is ignored. 

Suppose we periodically and adiabatically vary three gate voltages V±, V 2 , and V3. 
That is, we choose independent external parameters X 1 , X 2 and X 3 as the gate voltages 
Vi, V 2 , and V 3 , respectively, which allows us to control the scatterer in different ways. 
For instance, we can choose to adiabatically change Vi and V 2 with V 3 kept constant, 
i.e., Vi = V? + AVi smut, V 2 = V 2 ° + AV 2 sin cut, V 3 = V 3 °, (AVi, 2 < V\° 2 ), with lu being 
the slow frequency characterizing the adiabaticity. 

The allowed discrete fc's may be solved numerically from the dispersion equation, 
which results in time-dependent discrete energy levels when the dot undergoes an 
adiabatic cycle. Suppose there are totally N e = 2M electrons in the ring. Then at 
low temperature, they occupy the lowest M level, with two electrons in each level. The 
total current J takes J = 2£^ 1 / (0 = 1/(2ttt) E^i(0fe - with the subscript I 
labels the lowest M energy levels. In Fig. 3(a), we plot the total current, X7rrdefined 
positive clockwise, for 58 electrons as a function of the parameter AV 2 . The current 
is most sensitive to this parameter. We also plot the contribution to the current from 
the time dependence of k. This contribution is quite small. In Fig. 3(b) we keep the 
number of electrons constant at 10 and change the size of the ring. When the size of 
the ring is large compared to the size of the QD (L=100 and 200 compared to a QD 
size of 18), then, when increasing the size the current increases only very little. The 
most interesting behaviour is shown in Fig. 3(c). Here we plot the total current as a 




Figure 3. (a) shows the total current as a function of the contour AV2/V2, and 
contribution to the total current from dk. As we can see the contribution from the 
time dependence of k is quite small. The ring size is 40 (if the size of the QD is 18 in 
the same arbitrary units), and the number of electrons in the ring is 58. 

(b) shows the total current for the first 10 electrons for different sizes of the ring as a 
function of the size of the contour AV^/V^ . The current increases as a function of the 
ring-size for a fixed number of electrons. 

(c) shows the total current as a function of the number of electrons in the ring, for 
different sizes of the ring. Oscillations in the total current is clearly seen. We used 
AV 2 = 0.1K, , and AVi = 0.1V?. The size of the QD is 18. 



function of the number of electrons for different ring sizes. Compared to the persistent 
current [17] in the ring, there is no large odd/even effect where the current changes sign 
with the parity. However, in this figure we can see that the current can change sign (thus 
direction) as a function of the filling, and the effect is larger for a smaller number of 
electrons. This can be detected by measuring the magnetic field produced by electrons 
traveling around the ring. 

In conclusion, we have studied adiabatic scattering in a ring geometry and 
established the connection between the circulating (pumping) current and the geometric 
phases acquired by the time-reversed scattering states based on a pure geometric 
argument. This may be viewed as a generalization of persistent currents in a mesoscopic 



Geometric phases in a ring 



9 



ring, in the sense that the Aharonov-Bohm effect is a special case of Berry's phases. We 
have also performed numerical calculations on a quantum dot embedded in a small 
metallic ring. 
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